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We investigate the physics of frustrated 3-leg Hubbard ladders in the band limit, when hopping 
across the ladder's rungs (tx) is of the same order as hopping along them (t) much greater than the 
onsite Coulomb repulsion (U). We show that this model exhibits a striking electron-hole asymmetry 
close to half-filling: the hole-doped system at low temperatures develops a Resonating Valence Bond 
(RVB)-like d-wave gap (pseudogap close to (7r,0)) coinciding with gapless nodal excitations (nodal 
liquid); in contrast, the electron-doped system is seen to develop a Mott gap at the nodes, whilst 
retaining a metallic character of its majority Fermi surface. At lower temperatures in the electron- 
doped case, d-wave superconducting correlations - here, coexisting with gapped nodal excitations 
- are already seen to arise. Upon further doping the hole-doped case, the RVB-like state yields to 
d-wave superconductivity. Such physics is reminiscent of that exhibited by the high temperature 
cuprate superconductors-notably electron-hole asymmetry as noted by Angle Resolved PhotoEmis- 
sion Spectroscopy (ARPES) and the resistivity exponents observed. This toy model also reinforces 
the importance of a more thorough experimental investigation of the known 3-leg ladder cuprate 
systems, and may have some bearing on low dimensional organic superconductors. 

PACS numbers: 71.10.Pm;71.30.+h;72.10.-d 

I. INTRODUCTION 

It is a commonly held belief that the physics of the high temperature cuprate superconductors results from hole 
or electron doping a Mott insulator. This perspective has been reinforced by recent measurementsi*^ of optical 
conductivity which clearly show a large (but decreasing) Mott feature existing at an energy scale of order 2 eV 
to substantial dopings. One puzzle theoretically has been to understand how this feature remains relevant to such 
high dopings. We therefore propose a simple quasi-one-dimensional (ID) toy model which appears to capture the 
basic features of this physics and qualitatively provides a suggestive answer to this question. In addition, this model 
seems to be capable of differentiating between electron and hole doped systems which may qualitatively provide 
an interpretation of results seen by ARPES and provides a natural framework within which one might realize the 
formation of preformed (d-wave) pairs as an explanation of the mysterious "pseudogap"—. Thus, we have a (quasi-lD) 
realization of the point of view put forward by Anderson in 1987'*, that hole doping an RVB state leads to a (quasi 
d-wave) superconducting state. 

Ladder systems have been subject to extensive theoretical and experimental studies. As succinctly summarized by 
Dagotto and Rice^, it had been established by 1996 that the spin-i Heisenberg model shows alternating spin-gap/no 
spin gap behavior as one adds one-dimensional chains together, generalizing the Haldane conjecture as one increases 
the number of chains. Such behavior had been observed experimentally by Azuma et al^ in two ladder analogues 
of the undoped cuprate systems: SrCu203 a two-leg ladder exhibiting a spin gap; and Sr2Cu305 a three-leg ladder 
exhibiting no spin gap. About the same time, Schulai found that as the number of legs increases towards infinity, 
the spin gap disappears, raising the spectre that one might be capable of addressing the 2D limit of the Hubbard 
model-known to be an antiferromagnetic Mott insulator. 

Superconductivity in 2-leg Hubbard ladder systems was first predicted in 1994 by Sigrist et d& using a large U 
technique. Further evidence at strong coupling followed from Tsvelik and Shelton^. Starting from the opposite (small 
U) limit, Balents and Fisher»i2. showed that it was possible to classify all the possible ground states for the two-leg 
Hubbard ladder, developing a Renormalization Group (RG) technique used in collaboration with Linii to investigate 
the phase diagram of the Hubbard model away from half-filling. In this first paper, it was shown that the two-leg 
ladder at small U could also support the one-dimensional analogue of a superconductor-which was seen to have a 
d-wave like character in an excellent introduction to the area by Fisher*^. Furthermore, experimentally the telephone 
compound Sro.4Ca13.6Cu24O41.84, which is thought to be a hghtly doped coupled 2-leg ladder system, appeared to 
show some evidence for superconductivity at a pressure of 3 GPa as noted by Uehara et aU^. Solutions of the half-filled 
2-leg^'^^ and 3-legi^*iSi Hubbard ladders at weak U followed, and were generalized to the N-leg case where signatures 
of an antiferromagnetic ground state were founctf . 

Our starting point is the half-filled 3-leg Hubbard ladder, which is the simplest model which can hope to capture 
the basic physics of the cuprates. We will first review the interesting physics exhibited by this model as a function 
of doping and ask why it is that some qualitative features of the cuprates (such as electron/hole asymmetry close 
to half-filling) seem to be absent from this minimal model. Investigating the chemistry of the copper-oxide plane. 
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we find (following several band theory papers by Andersen and otherg-'^-) that an effectively repulsive next-nearest 
neighbor hopping should be added to any realistic (1-band) description of the cuprate systems. This causes us to 
revisit the three-leg ladder, adding a frustrating hopping term which qualitatively changes the physics for electron- 
doped cuprate systems, although unfortunately it also destroys the Mott insulating state at half-filling in the weak 
U limit (as found previously by Kashima et al in 2D)^^. While the cuprate superconductors certainly belong to 
the 2D limit of ladder systems (ie. oo), it is encouraging to see that many of the basic features are captured 

within this simple 3-leg ladder approach, allowing us to speculate on the nature of the N-leg system. While there are 
certainly many disadvantages to working with a ID model to describe 2D physics, it is difficult in 2D to go beyond 
the renormalization group approach and to provide a real microscopic theory of the low-energy fixed points It has 
thus been the approach of many to treat 2D gauge theories for which Senthil and Fisher— have argued the necessity 
of the observation of visons to any theory including spin-charge separation which has not as yet been accomplished 
experimentalljiS^. Such topological excitations are thought not be present in this simple (effectively ID) model. 



11. THE NEED FOR FRUSTRATED HOPPING 



A. ARPES and the doped 3-leg ladder 



1. A band picture 



If one works in the band limit (see Fig. 1) with tj^ > 0.3 t ^ U, Ledermann et al^^ found that the six Fermi 
points led one to consider 21 coupled 1 loop RG equations, deriving from Cooper (forward and backward) scattering 
and umklapp processes. Upon integrating out the high energy modes, it was found that 8 couplings corresponding to 
an effective two band problem between bands 1 and 3 at half-filling, as considered in great detail by Lin et ali^ for 
the 2-leg ladder, scaled to strong coupling, reaching the same fixed ratio sufficiently before the 1 loop RG procedure 
broke down. Upon bosonization of bands 1 and 3, they recovered the D-Mott state exhibiting both a charge and spin 
gap and possessing preformed Cooper (hole) pairs above the gap. (The name D-Mott is taken to emphasize that the 
superconducting order parameter changes of sign between bands 1 and 3). From the RG equations, one could estimate 
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FIG. 1: (a) 3-chains in the band basis: tx ~ t is the vertical distance between bands, 4t the band-width. Notice that at 
half-filling one has 6 Fermi points, and that vi=V3 < V2, where Vi can be read from Eq. 5 with t'=0; (b) Gaps appear in band 

i at energy « te " ' , leading to a sequential gapping of the Fermi surface. Points are shown on the corresponding quasi- ID 
(purple) or 2D half-filled Fermi surface; (c) The D-Mott state in real space: an RVB-like ground-state across the rungs of the 
ladder (spin and charge gap at i-filling). Per rung we have: 1 singlet and 1 unpaired spin; (d) The resonating Luttinger liquid 
state. The two outer legs feature one spin per site, one per rung being unpaired. At T=0 a Mott gap forms, which forces the 
spins on bands 1 and 3 to align ferromagnetically, opening a charge gap. 



the energy at which this gap opened as being of the order Ai « te f \ Having gapped this part of the Fermi surface. 
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the renormalized Fermi surface was now found to consist of just two points, and integration of the remaining high 
energy modes allowed the second band to subsequently flow to strong coupling, hence opening a charge gap at an 
energy scale of the order A2 ~ te u'^ . Thus it was seen that the three-leg ladder exhibited evidence of a truncated 
Fermi surface due to the differing Fermi velocities of different bands. This is illustrated in Fig. 1, where one sees 
that for open boundary conditions, the six Fermi points actually lie on the 2D Fermi surface as first pointed out by 
Lin et alii, although the perpendicular component of the momentum is only defined up to a sign. The hierarchy 
of energy scales resulting from the velocity differences between bands as pictured in Fig. 1(a), means that over a 
certain energy scale one would expect to have nodal excitations coexisting with a gap-like feature away from the 
direction-perhaps analogous to the mysterious pseudogap, and playing the role of "preformed pairs". Evidence 
for a truncated Fermi surface behavior has also been put forward in the context of quasi- ID organic systems such as 
(TMTSF)2PFg2^424i25i2SiSi although some recent evidence argues in favor of macroscopic phase separations^. 

2. The effect of doping 

The effects of perturbing such a state as a function of electron or hole doping are presented pictorially in Fig. 2. 
Notice that such a model possesses particle-hole symmetry meaning that a small addition of electrons or holes to the 
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FIG. 2: The first quadrant of the Brillouin zone, the dashed white line representing the 2D half-filled Fermi surface, (left) Un- 
doped ground state: Mott insulator; (center) Weak electron(hole) (purple (yellow)) doping ground state: D-Mott gap -|- "nodal" 
liquid featuring unpaired electrons(holes) on the outer chains; (right) Large dopings: approximate d-wave superconductor. 

i-filled Fermi surface would initially be expected to expand or contract the Fermi surface where the smallest gap 
exists. Thus, as pointed out by Ledermann et al^-, for an infinitessimal doping (at small U), one would expect the 
D-Mott gap to remain, with a small electron or hole patch of spectral weight to appear near {^,^). Such a picture is 
qualitatively similar to that espoused by Geshkenbein et al2S advocating the idea of bosonic preformed pairs (with a 
small spectral weight) in the corners of the Fermi surface coexisting with nodal (real) electrons, and the nodal liquid 
state proposed by Balents, Fisher and NayakJ?£ which rather featured the nodal liquid in terms of electrically neutral 
quanta called nodons plus insulating Cooper pairs. For the three-leg ladder this picture has been verified numerically 
by White and Scalapinc-^ who, using density matrix RG (on a t-J model), were able to dope a three- leg ladder and 
observed that the holes only occured on the outer legs (corresponding to doping the resonating Luttinger liquid state 
as pictured above in Fig. 1 (d)). Upon further doping, one would expect to drive the umklapp interactions irrelevant, 
at which point the D-Mott state was seen to be unstable to critical correlations in the d-wave superconducting order 
parameter at low temperatures. At higher temperatures, one would be left with metallic behavior. Similar truncations 
of the Fermi surface were found at large U by Rice et al^^ for a 3-leg t-J model, (despite the concern that |;-filling 
umklapp terms might become relevant in the outer two bands as U increased). 

3. Comparing to experiment: ARPES 

If one is willing to assume that ARPES results represent the properties of the bulk superconducting copper-oxide 
planes, despite the fact that it is a surface probe, one is immediately struck by the difference between the doped Mott 
insulating ladder systems and the cuprates. In particular, experiments seem to point to a clear difference between 
electron and hole doping. Spectral weight on the hole-doped side (Fig. 3 (left)) seems to first appear (in the normal 
state) close to the nodal direction along (§,§) as one might expect from the three-leg ladder, although perhaps there 
is some curvature of the Fermi surface observed. In contrast, on the electron-doped side, spectral weight appears first 
close to the (7r,0) and (0,7r) directions, to be followed at higher dopings by additional weight along (f ,f ). While there 
has been some recent speculation about the effect that annealing has on electron-doped samples (before annealing 
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they apparently do not show evidence for superconductivity^), we will assume that this has been performed in the 
recent experiment presented here, and ask how can such a bizarre behavior arise from a doped Mott insulator? Are 
we missing something in our simple Hubbard model? 
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FIG. 3: (left) 5=0.10 normal state hole-doped Fermi surface as extracted by ARPES^?:, notice the appearance of spectral weight 
close to (-1,^); (center) 5=0.10 normal state electron-doped Fermi surface as extracted by ARPES2", notice the absence of 
spectral weight close to (|-, ^); (right) 5=0.14 normal state electron-doped Fermi surface as extracted by ARPES^, notice the 
spectral weight both close to (7r,0) and close to (f , f) . Qualitatively only hole-doped features are described by Fig. 2 (center) 
above. The symbols PG, AF, NFL, FL, and SC in the phase diagram above denote PseudoGap, Antiferromagnet, Non Fermi 
Liquid, Fermi Liquid, and Superconductivity respectively. The NFL realm is embodied by a linear behavior of the resistivity 
versus the temperature. Region I has been a topic of recent debate, in La2-a;Sra;Cui-j;04 it has been argued that a spin 
glass ground state occurs, whose signature is enhanced by the addition of Zn on the magnetic sitea^, while clean YBa2Cu30j; 
samples appear to show physics one might attribute to a nodal metal—. 



B. Cuprate chemistry 

Band structure approaches to the high temperature superconductors do an extremely good job of explaining the 
physics at high energies, yet miss the crucial low energy physics that is thought to be responsible for the key low 
energy properties close to the Fermi energy. A particular difficulty for these approaches is a full treatment of the 
on-site Coulomb repulsion, thought to be responsible for the (Mott) insulating behavior of the 1/2 filled systems. 
Nonetheless, many important ideas do arise from such calculations which are able to derive not only an effectively 
direct hopping between copper sites (as used by the majority of strong correlation techniques), but additional effects 
resulting from the hopping via filled oxygen sites allow the contribution of unfilled copper orbitals, leading to the 
derivation of additional diagonal and second neighbor hopping terms which qualitatively change the physics of strongly 
correlated systems. The introduction of such terms has proven useful in several recent attempts to explain the ARPES 
spectrum of both electron and hole-doped cuprate superconductors within this communitji^fliiMiiiMMIil^, so it is 
worth reviewing the physical origin of such terms, before including an effective diagonal hopping. For simplicity, 
inclusion of second neighbor hoppings will be left for future work. The following is then a summary of the basic ideas 
of Andersen and collaborators-'^^. 

As one sees from Fig. 4, keeping only cr-bonding means that the orbitals of direct interest to us would be 
Cu{3da;2_j^2 ,4s} and O {2px, 2py}. If the distance between copper atoms on the lattice is a and R = {na,ma,0) then 
a choice {|Cud >= J2n,fn l^'^x^-y^ > ^'^'^^ l^x >= E„,m«"^|02p, > e*'' (^+(t.o,o))^ } leads to the matrix element 
<Cnd\Ox —2tpd sin(^l^), where the sign comes from the over-lap of the phases. In such a way, it is straightforward 
(in the absence of U) to generate a 4 x4 matrix of the dominant contributors. To express the resulting physics in 
terms of the copper atoms (as assumed by the majority of strong coupling treatments), one can imagine integrating 
out the fiuctuations to the oxygen orbitals then further integrating out the high energy fiuctuations to the unfilled 
Cu 4s level to generate an effective hopping between neighboring Cu 3d2.2_y2 orbitals. Following Andersen et allS as 
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FIG. 4: In the copper-oxide plane, while we are used to thinking about a 1/2-filled set of spins on the copper sites and excess 
holes or electrons hopping between them, it is useful to remember that hopping proceeds from one Cu atom to another via O 
atoms. The important matrix elements therefore are those between the nominally filled 2 Px,(y) orbital and the unfilled second 
electron orbital on the copper atom 3dj.2_y2. Virtual processes may also allow hopping between the unfilled Cu 4s orbital and 
the filled 2 Px,{y)- Matrix elements are denoted t, only cr-bonding of orbitals will be considered (following Andersen), and the 
relative phases generate the signs of each matrix element. 



reviewed in Appendix A, one obtains 

where s(e) = ^'^"^(^2^'^^''^''^ denotes the relative importance of hoppings between the Cu 4s and O 2p levels, u = 
^{cos{akx) +cos{aky)), and v = ^(cos(afcj,) — cos(aA:2:)) arises from virtual fluctuations between the 4s and 3d levels. 
This leads to the dispersion relation 

e « H 7== ( 1 - d,{eF) ~ i(cos(afca;) + cos(afcj^) ) - ^{cos{aky) - cos(afc^))^ 

\/d[eF) V ^ / ^ 

X (l + r(cos(afca;) + cos(afc.y)) + r'^{cos{akx) + cos(afcj^))^ + ..)) , (2) 

where a is the lattice spacing, d{e) = ^'^ {21%^'^'''' ^^'^ ~ 2{i+l{eF)) ' ^■^ expanding this series one can simply 
obtain the relative strengths for t,t',etc. in terms of a two-dimensional tight-binding model among the copper atoms: 
H = < e > —2t{cos{akx) + cos{aky)) + At' cos{akx) cos{aky) — 2t" {cos{2akx) + cos{2aky)). The leading dependence 
of each of these terms yields (in units of t« |(1 -I- r^), t'~ j and t"« |. Note that these results imply that 
because of fluctuations to the Cu 4s orbital, one generates an effectively repulsive hopping contribution (countering 
the hopping probability from two hops via the d-orbitals alone) . The sign of this term comes from the phase-shift of 
the cross-term tpstpd as expressed through the parameter v, resulting directly from the symmetry of the underlying 
orbitals on the lattice as derived in Appendix A. Inspection tells us that j ~ '^j while jr ^ \- 

In addition to the in-plane effects of electron or hole-doping, it is interesting to consider the out-of-plane effects. 
One wonders for instance if the size difference between substituted atoms might have a significant effect. In particular, 
would it be possible for t' to vanish close to half-filling, and grow with the doping? Andersen et al^^ have argued that 
the qualitative description above does not change markedly, but that the energy of the non-interacting Cu 4s band is 
modified via a hybridization to the apical oxygen atoms, such that inclusion of the apical oxygen atoms, metal atoms 
and the Cu ZA^^^^i orbital renormalizes its energy changing the effective strength of t'. The derivation of this result 
has been reproduced in Appendix B. 

C. Frustrating the 3 leg ladder 
1. A minimal model 

The minimal model including a contribution from this virtual hopping process is to include only the first term in 
the series, that corresponding to next-nearest neighbor hopping, t'. We will leave the t" (a double hop) for future 
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investigators on the grounds that it is only half the size of t', and thus should provide a smaller correction to the 
physics. We then proceed with the model detailed in Fig. 5, which written in terms of operators is, 

t -f 

FIG. 5: Inclusion of a frustrated hopping term. Hopping along the ladder (t) is treated in the continuum limit, hopping 
perpendicular to the chain {tj_) is of the same order and the Hamiltonian is diagonalized subject to open boundary conditions. 
The new term (f) enters with the opposite sign, and can naturally be included into this formalism. Notice that it is an 
effectively repulsive next nearest neighbor hopping-the sign arising from the phase difference between overlaps with oxygen 
atoms due to the integration out of virtual fluctuations to the unoccupied Cu 4s orbital. 



tdl^{x)dia{x + 1) + t±dl^{x)di+icr{x) 

- t'{dl^{x)d^+i^{x + 1) + dl^^^{x)d,^{x + 1)) + h.c. + C/n,|(x)n,x(a;)), (3) 



where N = 3 denotes the number of legs of the ladder, t represents the hopping amplitude along the chain, t_L the 
hopping perpendicular to the chain, t' the frustrated next- nearest neighbor repulsive hopping contribution, and U the 
(weak) strength of the on-site Coulomb repulsion. As in the unfrustrated case, we proceed by diagonalizing the U=0 




FIG. 6: The dotted lines represents the band structure (t=tx) of the three leg ladder when t'=0. Solid bands correspond to a 
choice t'= 0.1 t. Horizontal straight lines denote different 3-band Fermi surfaces. The lowest line represents the doping at which 
2-band umklapp processes dominate the physics (slight hole doping) (see also Fig. 9 (left)), while the dashed line represents the 
doping at which 3-band umklapp scattering is relevant (slight electron doping-at half-filling for even ladders) (Fig. 9 (center)). 
The half-filled location of the 3 leg t'=0.1 t Fermi surface is not shown but lies between these two lines. The uppermost line 
now corresponds to the filling at which the single band umklapp becomes relevant (sufficient electron doping-at a chemical 
potential which formerly corresponded to half-filling of the unfrustrated 3 leg ladders) (Fig. 9 (right)). 

ei = —2t cos(fc) + At' cos(fc) cos{kyi) — 2t± cos{kyi), (4) 

where, following Balents and Fisher—, we have first taken the continuum limit in the x-direction, before diagonalizing 
the resulting matrix in the y-direction, and kyi = We see from Fig. 6 that the effect of adding this frustrated 

hopping term is to break the degeneracy between the Fermi velocities of bands 3 and 1 , slightly raising the former 
and lowering the latter. The Fermi velocity (in the x direction in which the continuum limit has been taken) can be 




N f 
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simply found at any filling level (after solving for the value of kpi) from the equation, 

VFz = ^ U=feFi= ^2i-4i'cos (jif-^Ti^^ sin{kFi). (5) 

Due to the frustration of the lattice, the chemical potential of the half-filled Fermi surface decreases slightly. We 
continue to use the ID filling condition = n to denote the filling fraction, n, of the A^-leg ladder system, 

despite noticing that this condition now disagrees with the 2D filling condition '"'•J'"^^ — ri. Note that this 

definition means that, since we have electrons with spin, half-filling corresponds to 1 electron per site on average. 

III. THE PHYSICS RESULTING FROM FRUSTRATED HOPPING 



A. Interactions 
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TABLE I: The three leg ladder couplings, see Fig. 7, and derivation of their bare values from the Hubbard model. Note that 
spin umklapp terms are not present in the bare Hubbard model, but are spontaneously generated under the RG transformation 
if they are not disallowed under symmetry because the current operators belong to a high symmetry group. Here in terms of 
the Pauli matrices, ct", the definition of the (U(l)) charge current operator is J^ij = '^liis'^hja, the (U(l)) charge umklapp 
current is lhij= J2b,b' '^'^"^"'^'"'('Ljs' , the (SU(2)) spin current operator is J^^^ = I Ss.s' ^'LaO'L' V''>i3 and the (SU(2)) spin 
umklapp current operator is If^ij — ^ ^, V''iis(eo"°)ss' V''ijs where e — —ia^ as defined previously by Ledermann et al.— - Left 
and right movers have been identified at low energy as ipio- = e''''''^'ipRia(x) + e~''''''^ipLirT{x) with x = ja; see section III C. 
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FIG. 7; 3-band and 2d representations of the interactions taken into account in our treatment of the 3 leg ladder. Dark blue 
denotes Cooper processes, dashed dark purple forward and backward scattering and red umklapp scattering. See Table III 
for a detailed description of the different contributions. Here t'=0 has been taken for simplicity. The purple line in the 2d 
representation is drawn to note that the 3 leg ladder results also lie on a quasi-lD Fermi surface, in contrast to results when 
the number of legs grows. 

All interactions between Fermi points which conserve momentum up to lattice wavevectors are considered, as shown 
in Table I. For a picture of these allowed couplings, see Fig. 7. These four-fermion interactions arise from consideration 
of the onsite Coulomb repulsion. At one-loop order, these represent the only allowable terms, and the frustration 
of the hopping has the effect of reducing this list further. Generally at weak coupling, umklapp scattering processes 
are only relevant over a short range of doping, and away from this doping an oscillatory prefactor enters to negate 
such terms. The special doping level at which the umklapp processes shown in Table I are relevant in the absence 



9 



of frustration occurs at zero doping (or very close to half- filling) . The bare values of couplings are extracted directly 
from the Hubbard model. By writing these interactions in terms of current operators, as demonstrated by Balents 
and Fisheaiii, one has a controlled manner of writing the Renormalization Group (RG) equations by respecting the 
natural SU(2) spin and U(l) charge of the problem. Notice that spin-umklapp terms can be generated, upon scaling, 
from charge umklapps despite their initial absence in the Hubbard model. In the presence of frustration, electron-hole 
doping asymmetry rears its head to allow different umklapp processes to dominate the physics at 3 different dopings. 
At weak U, umklapp processes are no longer relevant at half-filling (for odd-leg ladders) so that the frustration of the 
lattice destroys the Mott insulator at half-filling. In the large U limit, it is expected that the domain of applicability 
of each umklapp process will become quite large, such that the three cases described below should gradually meld 
into one another, so there may be a substantial range of Mott insulating character in this limit. 

As there is reasonable convergence to fixed ratios of couplings well before the couplings become of the order of the 
band width (see Fig. 8) we do not consider higher order terms (which should be irrelevant). 



^ C [j/C ,3 
B □ 4C'',3/C''|3 

O B> 4fjc\^ 

B □ 4u'',33/c'', 

O .-8u''..33/c'', 
o o-c^/c'',3 

^--*-C 33/C ,3 
- - f 13/C 13 

O u\,,Jc\, 

A A c''33/c'',3 



lV° 10"' 10"' 10"' 10"' 10° 

FIG. 8: The frustrated 2-leg ladder RG at t=tx, t' = 0.1 t and ln(Y)=-21.8. Cyi/it — t') does not change signs under the 
renormalization group scaling, so that it is possible to follow the evolution of the ratios of the couplings as it scales towards 
strong coupling. Note that at strong coupling — ^ 1.11 so that the x coordinate is essentially the distance from the strong 

coupling values. We see that the 1 loop RG equations scale towards fixed ratios at a reasonable distance from the strong 
coupling limit. 



B. RG flow analysis 

As illustrated in Fig. 9, for the 3-leg ladder, we now find three special fillings close to the half-filling condition. 
We would like to stress that for the 3-leg ladder at weak coupling, half-filling is not one of these three special cases 
where umklapp scattering is relevant such that, although Mott physics is present in our system, it does not lead to a 
Mott insulator at half-filling if t' remains large with respect to U. Very recently, Kusunose and Rice-- have been able 
to capture both AF and a progressive Fermi surface development similar to that seen in electron-doped ARPES by 
inclusion of t' to a Hubbard model at U/t « 8. While it would be interesting to investigate the ramifications of an 
N-leg ladder, as sketched in Appendix C-which perhaps could be relevant to a discussion of the differing gap energy 
scales seen by recent Scanning Tunnelling Microscope measurements on the underdoped side of the cuprate phase 
diagram by Davis et alS2i, we will content ourselves here with a description of the physics of these three cases. We 
see clearly that the electron-hole symmetry of the original 3-leg ladder (see Ref. 16 (t'=0, 5 = 0, Vi?i = ^rF7i) and 17 
(t'= 0, (5 7^ 0, vfi 7^ Vi^a) for the RG flows) has given way in our more realistic model in such a way as to preserve 
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(0,0) 



FIG. 9: A schematic diagram of the three regions of interest to this paper for the 3-leg ladder close to half-filling. The 
introduction of a frustrated hopping term t' results in a varied phase diagram where Mott physics has different manifestations 
as a function of hole or electron doping. Here we have set t=tx = 10 t', and see that as a function of doping S away from half- 
filling (5 = 0), we find that the gapping of different regions of the Fermi surface leads to qualitatively different physics, (left) 
Hole-doping: T> Ai{vfi,vf3) metallic, T< Ai{vfi,vf3) D-Mott (ID d-wave RVB) -|- doped Luttinger liquid. Further hole 
doping will create power law d-wave superconducting correlations at low temperatures {Ai{vfi,vf3) ~ Ai); (center) Close to 
half-filled: 3-band umklapp terms are relevant. Possibly an antiferromagnetic precursor (at half-filling 4-band umklapp terms 
will be relevant for the 4- leg ladder although they appear not to generate an antiferromagnetic coupling); (right) Electron- 
doping: T> A2: metal, Ai{vfi,vf3) <T< A^: Mott insulator -I- metal (spectral weight near (7r,0)), T < Ai{vfi,vf3) Mott 
insulator + approximate d-wave superconductor. (A2 =A2 with V2 «f2, Ai (ufi, ufs) < A2 ) 



the interesting "preformed pair" pseudogap scenario on the hole-doped side and spectral weight well away from the 
(■|,|-) direction on the electron doped side, in quaUtative agreement with the above-mentioned results for ARPES. 
These cases will now be discussed in more detail. Set tj^ = t, and t' = 0.1 t. 

1. kpi + kpz = tt: hole doped (5 = -0.021) 

In this case, 2-band umklapp processes between bands 1 and 3 are relevant. This has the effect of driving the system 
towards an insulating 2-band fixed point punctuated by a metallic patch along the nodal direction arising from band 
2 as strong coupling is approached. In real space this corresponds to unpaired holes propagating along the outer legs. 
Remembering that the unfrustrated 2-leg ladder exhibits S0(8) symmetrj«i^ at strong coupling, one sees (Fig. 10) 
that the introduction of frustration, only breaks the symmetry between the Cooper channels cli and C33 but the same 
cast of couplings is relevant. For the case t' = 0.1 t, we find: 

g~c1^~ 4/f3 « 4c^3 « 4m^33i « 8u5'i33 « ^^331 « -\{cl^ + c^a) with ^ « — . (6) 

2 C33 VF3 

In passing, it is interesting to remark that, in the special case t = tj^ thought to be relevant to the cuprate supercon- 
ductors, one finds that bands 1 and 3 are quarter and three-quarter filled respectively. Although at weak U this is 
not expected to provide a relevant quarter filling umklapp term as such terms arise only at third order in U, (while 
the tree- level half-filled umklapp relevant here is of first order in U), one might expect to see new physics arise in the 
large U limit. 

2. 2 kp2 + kpi + kp^ = 2 tt; slightly electron doped (5 = 0.013) 

This case is the closest relative of the half-filled frustrated even leg ladder (where 4-band umklapp scattering is 
relevant). Here, the only relevant umklapp scattering processes mix all three bands. At the fixed point this leads to 
a complicated hierarchy, 

9 ~ ^/l2 ~ ^'^12 ~ ^"1223 ~ ^/32 ~ ^"^32 ~ 4m5'223 ~ ^'^12 ~ ^'^221Z ~ ^fl2 ~ '^^32 ~ ^Cgj. (7) 

The large number of relevant operators at the fixed point poses a difficulty for bosonization, as one finds that Klein 
factors from the relevant umklapp terms do not commute with one another nor with the Cooper terms c^2 ^^'^ '^32- 
This case is not of primary interest to this paper, and will be dealt with in Appendix D. 
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FIG. 10: Coupling ratios versus frustration: the frustrated 2-leg ladder at strong coupling {j^^ = 0.2,ln(Y) = —21.8). As 
t'^ 0, we approach SO(8) symmetry. Slightly away from this special point, one still finds insulating behavior. Notice the 
high degree of symmetry remaining at the point considered in this paper, t'=0.1t. Many authors have argued that the high 
temperature superconductors lie between t'=0.1t and t'=0.35t. On the other hand, this phase diagram may also be of interest 
to magnetically frustrated systems, as for the two-leg ladder one has a simple sign change of t', which manifests itself in an 
interchange of c^'^ and c^'^ . In this context a much wider range of t' may be experimentally accessible (perhaps umt' —f —t 
might be of interest to the pyrochlore systems although the 1 loop RG breaks for t'> 0.76). 



3. = f ■• electron doped (5 = 0.063) 

In this case, only band 2 umklapp processes are relevant in the weak U limit. This leads to a (Mott) gapping of 
the charge modes of band two at resonably high temperatures, with 

4, = 2uP^. (8) 

This effectively decouples the second band from bands 1 and 3 and this is followed at much lower temperatures by a 
(d-wave) superconducting transition (with gapful nodal excitations) . Above this superconducting transition, bands 1 
and 3 do not exhibit a pseudogap close to (■n.O). In real space, the onset of this Mott gap appears to imply (see Fig. 
11) a ferromagnetic alignment between electrons in chains one and three, with some propensity to double filling by 
a resonating orbital reminiscent of the pi orbitals of benzene rings. In the low energy effective description for the 




(a) (b) (c) 

FIG. 11: The Mott condition ■\p\^{x)'\p2s{x) = 1 corresponds in real space to a linear combination of electrons over the two 
outer chains of the three leg ladder. Here chain 2 is shown as a dashed line. At any rung along the ladder one finds that the 
formation of a Mott gap in band 2 corresponds to one of the pictures: (a) one electron in each of chains 1 and 3, with spins 
aligned; or a resonant state with (b) 2 electrons in chain 1 fluctuating to (c) 2 electrons in chain 2. That is, wavefunctions of 
the form, X/s'-^is'-^L (1 + (d|^ +d3s))|0 >, where s has the opposite spin from s. At half-filling with t'=0, the ground state 
of the resonating Luttinger liquid state of Fig. 1 (d) corresponds to a = or case (a) here. 

physics of bands 1 and 3 we are left with only 4 Fermi points. The lower temperature fixed point is characterized by 
the approximate values, 

5 = c^3 « 1.9c?i « 8.3cg3 « -1.4c?3 « -U.TS/fg with ^ = c^g = -cjg and /fg « 0. (9) 

C33 WFS 4 
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FIG. 12: (left) Coupling ratios versus frustration; the frustrated 2-leg electron doped (S = 0.05) ladder at strong couplings 
{jE^ = 0.2, ln(Y )=-21.8). (right) Coupling ratios vs. doping: the frustrated (f = 0.1 t) 2-leg ladder at strong couplings 
(■j5f7 ~ 0.2, ln( y)=-21.8). In both cases g is taken to be the largest coupling. Notice that in both cases, the absence of 
umklapp scattering means that the highest symmetry point differs from the SO (8) fixed point in three couplings: f^'g decreased 

by a factor of 2, and — are of the same magnitude. This turns out to be an S0(6) fixed point, as first realized by Schula^, 
which occurs on the geometrically frustrated side (t'<0) of the former diagram, but on the hole-doped side {S ~ -0.05) of the 
latter, and corresponds to the condition v^i = vfs- 



These ratios of couplings allow us to bosonize the remaining two-band contributions to extract the physics which we 
will show to exhibit superconducting correlations which dominate the physics close to the strong coupling fixed point 
at low temperatures (and is well compatible with a Mott gap in band 2). For comparison with the 2-leg frustrated 
ladder in the absence of umklapp terms see Fig. 12. We see that as a function of doping or t', there exists a special 
point in the absence of umklapp scattering analogous to the zero doping, t'=0 state of Fig. 10. This fixed point 
corresponds to the condition when vpi = ^fs, and possesses an S0(6) symmetry first noticed by Schulz^. The values 
of the couplings at this point are: 

5 = = ^11 = 8^3 = -c?3 = -8/f3 = 8c?i = -4c?3 and f^, « 0. (10) 



C. Bosonization 



Before proceding, it is perhaps useful to provide a brief introduction for the uninitiated into the world of bosoniza- 
tion. At low energies, provided the Fermi surface does not cut any of the bands too close to their minima/maxima, it 
makes sense to linearize the energy spectrum about each of the Fermi points. As such, we can replace the full band 
operator with local band operators whose momentum is measured with respect to the Fermi wavevector. By doing a 
partial Fourier transform, one arrives at the replacement, ipi^ — e^'^''^^ipB.ia{x) + e~^''''^''"tpLi^{x) with x = ja, where a 
represents the lattice spacing between sites and j the number of sites in a given displacement x. To cast these spin-ful 
fermionic operators in terms of bosonic operators, it is consistent to choose ■0pi s = ^''^ e'V^(-^^'^Pi+'"^'^»^~'^^p»+''^''i)^ 

' ' v27ra 

where P = (+,-) = (R,L) and s = (+, -) = (up, down); 77,. s is a Klein factor taking into account the fermionic anti- 
commutation relations between different bands ({^yj.s, Vjs'} = '^^i.jSs,s')', and {(f>pi, Opi) are conjugate bosonic variables 
corresponding to the charge degrees of freedom (and defined as (t)pi = , and (^o-i, O^i) their counterpart in spin 

(defined as (ji^i = . That is, O^i = (j)Lui - (j^Rui and cjy^i = (l)Lut + (^Rvi satisfy tt/jj] = i5{x - y)5ap5ij, 

where tt^^ — dOf^i. Fermionic anticommutation relations within the band are taken into account by means of the 
commutation relations amongst the bosonic fields, [(j)+ai,<t'~i3j\ = ^"^"'^/"^ and [4>±ai{x), (l)±iij{y)] — . A 

more complete introduction has been provided by Shankas^. 
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1. The hole-doped 3-leg ladder 

We have already seen that band 2 remains metallic in this case (forming a doped Luttinger liquid), so the question 
of interpretation of the physics of the two outer bands is of import. This can be addressed at both small and large 

frustrated hoppings. 

a) Moderate t': For general values of t' we ean extract Luttinger liquid parameters, after diagonalizing the kinetic 
contribution to the interaction. A diagonalization of the charge degrees of freedom for a more general model (the 
electron-doped case) has been treated in Appendix E. The current case corresponds to a choice g=0 in Eq. 51, which 
has the effect of replacing Vd± in Eq. 52. The Luttinger liquid parameters are then given by Eq. 24-26 subject 

to the choice g = 0. There are two contributions to the fixed point Hamiltonian density, non-umklapp terms, 

Hl,t = j^^Wii cos(%/2/30<,J + c?3Cos(x/2/3<^<,3) - 4c?3 cos(x/2^(e^, - epMcos{V2^ie,, - 6,,)) 



+ cos(V27r((/)^, - 0^3))} - f^3(cos(V27r(0p, - 6lp3)){- cos(V27r(6'^, - 0^,)) 

+ cos(V2^(0,, -(/.„3))} + 2cos(V2^(6'p, - 6*^3)) cos(v^(0<., + (/'<.3)))}, (11) 



which simplify due to the symmetry 0^3 — Ac^^ to yield. 



2 rkr 

I 

+2|c?3|cos(v^(6'p, -6'p3))cos(V2^0^Jcos(V2^(/)^3)}, (12) 



"^"^ = (2^^ 2 ~ + ^(2cos2(x/2^</,.3) - 1) 



and umklapp terms, 

H^ot = (^;^{~"i33i cos(\/2^(.^pi + (l)p3)) (2cos(\/2^(.^<,i + </><,3)) + cos(\/2^(</><,i - (f)a3)) + cos{V2^{e„i - 6^3))) 

+471^33^ cos(\/27r(0pi + (j)p3)) (^cos(\/27r(6'CTi - 6*^3)) - cos(V27r(0<Ti - '/'(ts))) 

-16mi;^33 cos(A/27r(6'pi - 6*^3)) cos(A/27r(0pi + (/)p3))}, (13) 
which simplify due to the symmetry U1331 = 4ui33i to yield, 

H^L = ■^^^{2|<33i|cos(\/2;^(<^pi +0p3))cos(\/2;^(</.^i))cos(\/2;^(<^<,3)) 

+8K133I COs(V2^(0pi - ep3)) C0s(V2^((/>pi + <j)p3))}. (14) 

To minimize H^^^ + H^^^, the ground state then pins 

\Phx4>a^ = 7rn, \/27r<?!)<73 = ttto, \/27r(^pi - ^pa) = 7rp and \Plv:{<^p\ + ^ps) = 7rg, (15) 

where {{n, rn,p}, g} j- are integers whose sum is even. The first and second conditions are already expressed in 
terms of eigenvectors, the third pinning condition finds expression in terms of the diagonal basis as Eq. 27 with the 
substitution V(j+ — > v^, while the fourth finds expression as. 



^1 + («, + yrr^ ^\ + {v, vtt^ 



■p_b I = ■Kq. (16) 



Together, these two conditions have the efi'ect of pinning the charge degrees of freedom of the outer two bands such 
that it forms an insulator (remember that band 2 corresponding to nodal excitations remains metallic), and the two 
spin conditions mean that one additionally has a spin gap so long as the couplings in Eq. 12 and 14 remain finite 
(as they do for the half-filled 2-leg ladder for all t' we can access as shown in Fig. 10). This means that similar to 
t'— > 0, even though S0(8) symmetry is broken, quasi-long range superconducting order is not possible for this case, 
and such correlation functions will decay exponentially, leaving only 2(kpi +kp3) critical charge density wave (CDW) 
correlations. Similar to t'^0, this hole-doped state therefore possesses a d-RVB like pseudogap feature punctuated 
by a nodal liquid. 
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b) A special limit, t' 0: When t' remains small (see Fig. 10), the difference between the couplings and C33 is 
small, so it seems reasonable to work in terms of the SO (8) basis, but introduce different couplings, g, g+ = -C33 and 
g_ = -cli to annotate the effect of perturbing this symmetry with t'. Then focusing on bands 1 and 3, one finds, 

H — Ho + Hkin + Hpot, (17) 

where 

^0= J2 ^iid^'^.y + M')^ (18) 

i—1.3,iy—p.(T 

the contribution of the interactions to the kinetic part of the Hamiltonian density is, 

Hk.n = ^ Iff- {{dcj^.J' - (vr.J^) + .9+ ((90.3)' - (7^.3)')) + 2g((<9</.pJ(a0p3) ~ 7rp,7r,3)} , (19) 

and since for small t' we have approximately (5+ + .9-) = 2g, if we re- write this in terms of bonding and anti-bonding 
operators between the bands (0i/± = '^{'Pi^i i (pi^a), {v — {p,cr))) and make the assignments, {(pi, 9i) = (0cr+, ^'ct+), 
i(f>2, 62) = {(t>a-, ^'o-_), (03, O3) = {(t>p^, 9p^), and (04,04) = (0p_ , Op_), this simplifies to, 

4 

HMn « ^ E - ^l) + (501902 - ^1^2) . (20) 

In this basis, the pinning potential term becomes: 

^pot-J^j^i^ E cos(/30a)cos(/30;,)j +^^__pisin(/30i)sin(/302). (21) 

It must also be noted that one additionally generates a term, 

5Ho = ^(9037r4 -I- 7r3904 -t- 90i902 + 7ri7r2). (22) 

As g+-g- oc (yf3 - "vpi), the second terms of Eq. 20 and Eq. 21 are approximately proportional to t', so that one 
retains an approximate SO (8) symmetry: 

^ = E(^«^'^-'^-) - f (E V-lrVa)' + 0{vpi - VF,), (23) 

a a 

where ip\ = (V'ija' V'la)' runs from 1 to 4, and is a Pauli matrix, as defined previously by Lin et at 

2. The electron doped 3-leg ladder 
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(a) General t': First we need to re-express the fixed point Hamiltonian in terms of canonical fields. This procedure 
commences with diagonalizing the kinetic contribution to the Hamiltonian and extracting the corresponding Luttinger 
liquid parameters. Then we can re-write the interacting part of the Hamiltonian in terms of this diagonal basis, which 
allows us to see the low energy pinning of the (classical) ground state of the system. Armed with this information 
we can then compute the correlation functions for a number of order parameters-here we demonstrate this on the 
superconducting order parameter which is found (see Fig. 13) to be dominant at low energies, and find a power law 
decay with a non-universal exponent in the region of phase space applicable to this method. For a derivation of these 
results, the reader is referred to Appendix E. The Luttinger liquid parameters are given by. 



\ Z TT V TT Z TT 



'..1 + ^9_^^ ^ '-rC-^^^r + (^)^ and u., = 1,.^ ^(^] (24) 



2 TT \ TT 2 TT y \ Att 
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Kp±b 



VFI+VF3 I 
2 ' 



i-f)±V(i-f)^( 



7 + {^7 



and Kr, 



(25) 



where we have extracted the proportionaUty of c^j 
density is given by, 



-g Ypi for i=l,3. The total kinetic contribution to the Hamiltonian 



kill 



Wpib) + 



P"^" 1=1,3 



U„iKn 



2K„ 



(26) 



The remaining contribution to the fixed point Hamiltonian density is given by the non-umklapp terms H^g^ as 



detailed in Eq. 12 above. To minimize this energy, the ground state then pins \j2'n^a-^ 



2TT(j)a 



and 



^/27^(^^PJ — 9p^ ) = np where n,m,p are integers whose sum is even. The first two conditions are already in terms of 
the basis of eigenstates. The latter finds representation in our transformed basis ((tt =f g) "-^^7/^^ — Vd^^)- 
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/27r 



"d+ 



Vd+ 



l + {vd+ + Jl + vj,) 



7p-b 



"d+ 



[Vd^ 



np. 



(27) 



In order to compute correlation functions, we must express in terms of our band picture the types of order parameter 
likely to be relevant. One candidate is certainly d-wave superconductivity, as previous studiesiSiiS on a 2-leg ladder 
system have shown this to dominate away from the half-filling umklapp processes. If we then imagine forming a 
singlet between chains one and two concurrently with a single between chains two and three (see Fig. 1 (c)) and 
re-express in terms of the band picture, we see that such an operator does not depend on band two: 



A|o > = ^(^T^i - 4i4] + 4Ai - 4ii4] + 'i-cOio > 

Using Appendix E, this order parameter is then found to decay as. 



<Ai(.t)A3(0)>|--^(^) ' V ' V . (29) 

(naV \ \x\ J 



Since the unfrustratcd two- leg ladder result is known, it is useful at this point to imagine setting vpi = vps (not too 
far from half-filling). In this limit, we recover (^a)^ (j~j) ''^ with 



K 



p+b 



V-^ TT^ 2 271- 
\ y'^ ^ IT ) 2 2ir 



(30) 



This has to be compared with the d-wave superconducting correlation function cx r (^'^) where = ^g^~^ ^ ~ 

Kp+b^"^ obtained by Schula^ for two coupled Luttinger liquids. 

For this fixed point, we found that the bare values Mpi = 0.1096, = 0.2017 and stopping the RG flows close to 
strong coupling gave g = 0.012 and f^g = 1.7x10"^ so that Vd+ = -83.3 was non-universal and at this point, 

o / 1 \ 0.253Kp+b-l-0.247Kp„b „ / 1 \ 0-523 

|<Ai(x)A3(0)>| = ^ - = . (31) 

(7ra)2 \\x\J (7ra)2 \\x\ J 

Alternatively, one can form a charge density order parameter as (for the 2-leg ladder beginning in a chain picture): 
OcDW = ni - na = ^{d\^di„ - dl^d^^) = ^{ti}\^{x)ijz^{x) + %pl„{x)iji„{x)), (32) 
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whose (leading) correlation function < {Ocdw{x))'^{Oc'dw{0))'^ > contains a term proportional to 



cos(2(fci^i + kF3)x) ( y 



-(1+ 



0+1 



(33) 



It is of interest to note that while in the electron-doped case, kpi + ^f3 > ^, ioi the earlier considered hole-doped 
case where this correlation function can become relevant, ki?i -I- kp^ — ^, permitting a physical description of the 2 
(kfi -|- kiT-a) oscillation as being periodic with the lattice-that is repeating every lattice spacing. The analogous 3- leg 
ladder CDW definition of Ocdw in real space corresponds to the linear combination, 



CDW — 



E 



nicr + nsa 



n2a 



(34) 





3.25 




3 




2.75 




2.5 


c 

OJ 

c 


2.25 


o 




CI. 
X 


2 


OJ 






1.75 


o 






1.5 


a 






1.25 




1 




0.75 




0.5 



-^d-wave s.c. 
-Ocdw 



-0.75 -0.5 -0.25 





t'/t 



0.25 0.5 



0.75 



2.75 - 

^ 2.5 - 

jf 2-25 - 

I 2 - 

I 1.75 - 

I 1.5 - 

d 

^ 1.25 - 

8 1 - 

0.75 - 



0.5 ^ 
-0.30 



band loss 



0-- 



-0 O'- 



-Od-wave s.c. 
Ocdw 



band loss 



o ' cm ' 



-0.20 -0.10 0.00 0.10 

6 doping away from 1/2 filling 



0.20 



FIG. 13: Correlation functions (a superconducting order parameter with d-wave symmetry, and a charge density wave order 
parameter) as a function of: (left) t' at an electron doping of (5 = 0.05; and (right) doping at t'=0.1 t. Shown is the two-leg 
ladder result, < O{x)O{0) > ~ j^^, as derived in the absence of umklapp scattering. 

(b) Another special point: Using the basis introduced in section CI. , it is not difhcult to show that the SO (6) 
fixed point of Eq. 10 is identical to the S0(8) fixed point, save for the depinning of the symmetric charge mode (j)p+, 
which is the reason that quasi- long range superconducting order is stable about this point. At this point one recovers 
Eq. 23 with a now running from 2 to 4. 

IV. CONCLUSIONS 

In this paper we have investigated the phase diagram close to half filling of the frustrated 3-leg Hubbard ladder in 
the band limit. Within this weak coupling approach Mott physics is found to be relevant at three special dopings. 
In the normal state we have discovered physics quite reminiscent of that observed experimentally using ARPES. 
Specifically, the addition of a frustrated hopping, motivated by the actual chemical structure of the cuprate planes 
has the effect of breaking electron-hole symmetry to (once more) promote a d-wave RVB state with gapless nodal 
quasiparticles on the hole-doped side; and a Mott insulating (■!,§) feature coexisting with metallic physics at the anti- 
nodal points, seceding to d-wave superconductivity at low temperatures on the electron- doped side. Extensions of this 
toy model to an N-leg ladder system, as outlined briefly in Appendix C, would suggest that on the underdoped side 
(of the hole-doped cuprates) the d-RVB state may coexist with a partial d-wave like superconductivity, originating 
from different regions of the Fermi surface. As the pseudogap (d-RVB state) disappears, one would expect that the 
quasiparticles begin to cover the entire Fermi surface, so that only after the complete disappearance of the pseudogap 
might one recover a truly normal metallic Fermi surface. Unfortunately, we have not been able to show that such 
a model supports an antiferromagnetic ground state close to half-filling, unlike that seen when t'=0 in the large N 
limit. This curious feature is likely an artifact of the consideration of U < t' here, although it could also arise should 
the value of t' change greatly as a function of doping. 
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In closing, it is perhaps useful to additionally stress the link of these unusual ladder systems to the plethora of 
interesting physical systems now arising in the growing area of frustrated materials. A particularly common structure 
arising in magnetic systems is the pyrochlore or corner-shared tetrahedral lattice as pictured in Fig. 14 (a). It is 
well-known that such a lattice finds a simpler partner in a mapping to the checkerboard lattice (Fig. 14 (b)). The 
frustrated ladders with which we have been playing (if one reverses the sign of the diagonal hopping) are equivalent to 
a modified checkerboard lattice (Fig. 14 (d)), so one might hope as this hopping becomes large to recapture some of 
the interesting physics emerging from these systems. As we can see from Fig. 10 and 12, we can access a much larger 
range oft' than could correspond to the physical case of the cuprates, but cannot reach the limit t— >-t' whereupon the 
upper band-width would become zero. Indeed, as we approach this limit, the 1-loop RG scheme breaks down before we 
are able to see the metal-insulator transition expected from 2D studies by Kashima et al^^:, and one is left wondering 
if the superconducting correlations survive as the dominant contribution in this regime. Were this the case, one might 
be able to resolve whether the superconductivity of LiTi204, which possesses a Tc of 13.7 K^S., has a conventional 
electron-phonon origin, (perhaps resulting from lattice distortions along the [111] direction in analogy with the charge 
ordering mechanism of A1V204^), or is more unconventional in nature. It should be possible to access this regime 
within the auspices of currently available numerical RG techniques (such as density matrix renormalization group) , 
and remains an open and interesting extension of this work. 




(b) (d) 




FIG. 14: If one considers t' to be an additional hopping parameter with the same sign as t (ie. t' = - a t, a > in Eq. 4), 
then one should be describing a frustrated hopping term similar to that found on a geometrically frustrated lattice. There have 
recently been several experimental systems found demonstrating interesting physical properties due to the large degeneracy of 
states on a pyrochlore lattice. While the lattice we study here is more complicated (and therefore less frustrated) than the 
pyrochlore lattice, there are some similarities, and one might expect to see similar physics resulting in the limit close to t— > t'. 
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APPENDIX A: REDUCTION OF A 4-BAND MODEL TO PRODUCE t' 

Following Andersen et ali^, we can write from Fig. 4 the 4x4 matrix shown in Table II. Making use of the Lowdin 
procedure to integrate/downfold the Hamiltonian, 

= Hw - [H" - ^]JfHjH' , (35) 

one reduces the 4-band model to the effective two-band model shown in Table III. As an example of this procedure, 
we compute 

= _(_2tp,sin(^))^(-2t,,sin(^)) - (2t,, sin(^))^(-2tp, sin(^)) 



^tpdtps f cos{aky) — cos{akx) 



(36) 
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\Cud > 


\Cu4s > 




o« 


< Cud\ 


Ed 





-2Vsm(^i|^) 


2Vsm(^) 


< Cu4,s 1 







-2tp, sm(^) 


— 2tps sm(^) 


< 




-2tps sin(^) 







<Oy\ 


2tpa sm(^) 


-2tps sm(^) 








TABLE 11: The dominant contributors to conduction between copper atoms according to band structure. 



where we see that the phase shift of the overlap between the d-orbital and the p-orbital manifests itself in a sign 
difference relative to the phase of overlap of s- and p-orbitals. This symmetry property of the orbital overlap is 
responsible for the change of sign of t' relative to t in the expansion below as diagonal terms are created from the 
expansion of the square of this term. Further downfolding of the high energy 4s level then yields an effective one- 





\Cud > 


\CU4,s > 


< Cud\ 













TABLE III: The effective 2-band model. Here u = ^{cos{akx) + cos(afcj;)) and v = i(cos(afcy) — cos(afca;)). 



band description which may then be compared with that used in most strongly correlated electron approaches. One 
obtains 

H = ed + — - — 1 - u - — - , (37) 

e-Cp V l~u + s(e)J 

where s(e) — (2f'"'^)2^''"' denotes the relative importance of hoppings between the Cu 4s and O 2p levels. 

In order to head towards a t,t',.. model, it is useful to examine the area close to e ~ ep. Solving the Hamiltonian 
for the energy is equivalentiSi to solving = —d{e) + 1 — u — where d{e) — ^'^^^2t^dyi'^''^ ■ Then expanding 

to linear order in d(e) and assuming that the energy dependence of s(e) is weak about the Fermi energy (if ep lies 
symmetrically between ep and eg then s(eF) — 0), it is not hard to express the energy as, 

e^ep + {diep))-' (^1 - diep) ~u~ j^"^) ' (^8) 

where r = 2{i+l(ep)) becomes part of a natural expansion parameter to yield approximate square lattice effective 

hopping parameters. (Note that diep) = — t — ~ xhen, 

e « + j^f , ( 1 - d{eF) - \{cos{ak^) + cos(afcj^) ) - ^{cos{aky) - cos(afc^))^ 
Va(e-F) V ^ / ^ 

X (1 + r(cos(afca;) + cos(afcy)) + r^(cos(afc2:) + cos(afc2^))^ + ..)). (39) 



APPENDIX B: DERIVATION OF OUT-OF-PLANE CONTRIBUTIONS 



If one accepts that the energy of the non-interacting 4s Cu band is modified by hybridization with the apical oxygen 
atoms-in turn affected by the metal atoms and the 'iA^y2_i Cu orbital, one generates a picture as seen in Fig. 15. To 
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estimate the new energy of the apical oxygen, it might be necessary to consider this hybridization with neighboring 
overlapping atoms. To calculate such a correction due to physics outside the Cu-0 plane, following Andersen et aU^, 
we start with a 6-band picture as shown in Table IV. Proceeding with Lowdin downfolding of the 5 contributors 



M 




M 



FIG. 15: The out of plane contributing orbitals. 





\Cu4s > 




\02p, > 


\M > 


\02p^ > 


02py 


< Cu4s 1 







-2tscsm{2^) 





-2tps sin(^) 


-2tps sin(^) 









^2t,,2 sm{^) 





2t^^2 sin(^) 


-2^2 sin(^) 


< 02p, 


-2t,csin(^) 


-2t,,2sm{^) 












< M\ 








-tcAie 2 










< 02p^ 1 


-2tps sin(ii|^) 


2tp,2sm{^) 













< 02py 1 


^2tps sin(^) 


2^2 sin(^) 











ep 



TABLE IV: The out-of-plane bands. 



yields. 



(40) 



The first term is easy to calculate, H34{H^^ — e)^^ H43 — while the others require the inversion of a 4x4 matrix 

and some approximations. One finds, 

(2i,,2 sin(^))2((e, - e){ep - e) - {2t,sr{sm\^) + sm\^))) 



{es - e){ep - e)(e,2 - e) - (sm'{^) + sin' {^)){{2t 2 - e) + (e, - e)(2V2)2) 



, (41) 



and 



(2t,, sin(^))^((e,2 - e)(6p ~ e) ~ (2V2)^(sin^(^) + sin^(^))) 
(e. - e){ep - e)(e,2 ~ e) - (sin2(2|.) + sin2(^))((2tp,)2(e,2 - e) + (e, - e)(2V2)2)' 



-4te,2t,,(sin(^))2((4tp,V2)(sin^(i^) + sin^(^))) 

(e, - 6)(ep - 6)(6,2 - e) - (sm'i^) + s\n'{^M2t^,)He,2 - e) + (e, - 6)(2V)2) 



(42) 



(43) 



20 



Then, if one assumes that t^^2{^s — cf) << tlpi^F — Cz^), it is justifiable to drop the second half of the last term of 
the denominator. If one then averages over ]ixXyXz in hybridizing the bands such that < sin'^(^) i, these four 
terms add to: 

{{2t,,2f{es -€)- {2t,,f{€ - £,2)) - {{2tp,f{2t,,2f + m,,2tsctpstp,2 + (2t,e)'(2V)"} 
2((e, - e)(e - e,2)(e - e^) + {2tp,)^{t - £,2)) 

_ -{{2tps){2t,,2) + {2tsc){2tp,2)f 



2(e - e^2)((e, - e)(e - ep) + (2ip,)2) 
l + =^l'(2i,,2)2 



)2 1 



(e.-e)(e 



(44) 



where the first two terms in the numerator have been assumed to approximately cancel reflecting an approximate 
equality of the likelihood of hopping from the Cu to the apical O site or vice versa for the 4s and 2>d^z2_i orbitals 
respectively. The latter term in the denominator can then be defined as to reproduce the renormalized energy of 
the Cu 4s electron where the factor of 2 comes from there being 2 "apical" oxygen atoms, 



e,5 



2{tsc? 

ep - f-c 



(45) 



where = Cr 



pi2 



This dependence on out-of-plane physics makes many predictions in 



terms of how one might increase the maximal Tc of a sample if it is related to the strength of t', as noted in the 
interesting paper by Pavarini et aU^, although one might like to take such claims as "the fact that Tcmax drops from 
92 K to 50 K when Y is replaced by the larger cation La in YBa2Cu307" with a grain of salt in light of results by 
Wada et al^. 



APPENDIX C: A LARGE N PSEUDOGAP PICTURE 

If we consider increasing the number of legs of the ladder system, what expectations might arise as to the nature of 
the pseudogap from this picture? One might expect a shifting spin-liquid contribution to the gap (shown in Fig. 16 as a 





angle-dependent d-RVB 
+ nodal qp's 




BCS-like 



weak s.c. 



nodal Mottness 




d-wave 



FIG. 16: A pseudogap with k-space angle-dependent physics would be expected to arise for ladders with both frustrated 
hopping and a large number of legs. In red we see the effects of Mott physics as a function of doping which define the type of 
quasiparticles allowed as a function of doping. One sees that the lower temperature establishment of a d-wave superconducting 
gap would only occur in regions of the Fermi surface which have not already developed a gap due to umklapp scattering. Thus 
on the underdoped side of the phase diagram shown, one would expect to find both a pseudogap and a superconducting gap at 
different position of the Fermi surface. Sufficient hole doping (left) allows the Mott physics to contribute at the maximal gap 
value corresponding to the superconductor just before the Mott physics would vanish entirely, presumably leaving a normal 
d-wave superconducting gap over all the Fermi surface, and metallic behavior above. On the electron-doped side (right) one 
would expect to see a pseudogap vanishing abruptly as a function of doping and of substantially different character. 
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large red gap), coexisting at low temperatures with regions of d-wave superconductivity (shown in Fig. 16 in blue). At 
high doping levels (above those at which umklapp scattering was relevant) one would hope that interference between 
holes would destroy the coherence of the superconducting order parameter. Within such a picture, it is interesting to 
note the presence/absence of nodal quasiparticles on the hole/electron doped side which would naturally account for 
at least the low temperature difference of the resistivity exponents. 



APPENDIX D: BOSONIZATION OF 3-BAND AND 4-BAND UMKLAPP TERMS 

1. In the absence of frustration 

As mentioned in the introduction, the half-filled Hubbard ladder has been found to exhibit no spin gap for the 
3-leg ladder, and 3-band umklapp processes are not relevant. The 4-leg ladder possesses different physics: a spin gap 
at half-filling, while 4-band umklapp processes are not relevant. However, in the limit of large N (N even), it has 
been showni^ that the high degeneracy of umklapp processes about half-filling has the effect of driving all 4-band 
couplings, c'', j-^, c'^,T^, vl^-,t^ and u'^, r^, relevant. These terms undergo a large renormalization as, for any fixed i, 

^ ^ ' jkkj^ jkkj^ jkkj jkkj' o o 7 .y ; 

the condition kpj + kp^ -I- kp^i + kp^ = 2-k has N solutions, and enters the RG equations in such a way that each 
contributes. Upon bosonization of these surviving terms, Ledermannii found them to describe an antiferromagnetic 
state (with no spin gap). At half-filling, in the presence of frustration, we hoped to recover this antiferromagnetic 
Mott insulating state, but as we shall see below, it does not seem to arise. Indeed, the large N limit of the half-filled 
even leg frustrated Hubbard ladder likely has physics similar to that shown by the 4-leg ladder, while interpretation 
of the odd-leg ladder results is clouded by the non-commutativity of the Klein factors. 

2. In the presence of frustration 

We first attempt to bosonize (part a)) the 3-leg ladder for the case described in HI B 2, close to half- filling when 
its 3-band umklapp processes are relevant. While much of the fixed point can be bosonized, these 3-band umklapp 
processes do not allow us to proceed further and determine the nature of the physics. This is not the case for even-leg 
ladder systems at half-filling, as we show (in part b)) by bosonizing the counterpoint of 3-band umklapp terms for 
this system-4-band umklapps. However, it is not clear that the physics of the frustrated 4-leg ladder at half-filling 
should yield the same physics as seen at the special 3-band filling of the 3-leg ladder, as shown in Table V. Inspection 
of the RG flow of the frustrated 4-leg ladder equations further shows us that the 4-band Cooper scattering terms 
remain small, so that we do not recover the antiferromagnetic Mott insulator for the frustrated 4-leg ladder. But 
what happens if we go further, to the large N limit, where the unfrustrated ladder showed this interesting physics? 
As shown in part c), it is unlikely that in this limit the 4-band Cooper scattering processes would grow, as one no 
longer has a large degeneracy of band pairs lying along one's Fermi surface. 

a) N^3. If we first leave aside the 3-band couplings, the close to half- filled case can be bosonized as 

^ ^ X! lfi2i9<l>pid(l>p2 - T^piT^p2) + -^{d(l)aid(j>a2 ~ TTct j 7r<T2 ) | , (46) 

and 



TT 

1=1.3 



Hpot = J2 { cos{V2^{ep, - ep2))((-44 + <2) cos(^/2^(^^,, - 9,,)) - (4^ + c^,) cos(\/2^(0,, - ^,2))) 

—20^2 cos{V27r{6pi — 9p2)) cos(\/27r(0 

at \ c^a2)) + 2/^2 cos(V2^(0,, - e„2))cos{V2^{c^,, + 0,2))}, (47) 
which, upon choosing (2Tra)^g — — '^32 ~ ~'^'^32 ~ /12 ~ '^12 ~ ~'^^i2 simplifies to 

Hpot « 25 ^ {cos(\/2^(6lp, - 6ip2)){cos(\/2^(6l<,, - 61,2)) - cos(\/2^((/),,: + 0,2))} 



i=l,3 



-Hcos(%/2^(6i<,, - 9a2))cos{\f2^{<t>a, + 0^2))}- (48) 

Returning to the 3-band terms, one recognizes that the Baker-Hausdorf commutation relation vanishes for 
all these fields with the exception of the expansions of terms with two operators from the same band pos- 
sessing the same spin such as V'ij2s''/'i2s- Bosonizing this one finds a term « g-jV47r0R2ag-»\/47r0i,2s — 

g-zv^0R2»-iv^0i2ag|[-i>/47r0R2<,,-iV4ir0_L2s] _ g-iV4ir0_R2s-jV47r0i2s g^j^ _ _jg-iV47r0H2s-iv'^0l,2s 
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The umklapp scattering terms are thus bosonized as, 

J2mi)^ ^ {»72s'7is'72s'73s{«223 + M223) cos{(j)p++ + s(t>„^) cosiOp- + sd^ — ) 

-8^2213 cos(6'p++ + s6'ct_ ) cos((/!)p++ + - ?7isry3s{2 cos(^p++ + s0cr++) sin(6'p_ + 56*^-) 

+«223 - M223) cos((/)p++ + s^<,_) sin(6»p_ + s6»<,_)}}, (49) 

where we have defined cj)p++ = y|(2^p2 + <^pi + <i>p3), 0cr- = ~ ^<^3), Op- = ~ ^ps), 0^ — = 

-\/f^(2^<T2 — ^ct3 — ^cri), and ?iio-++ = y^(2^o-2 + i'ai + ^0-3 )• And a nasty surprise awaits. In order to be able 
to simultaneously diagonalize the Hamiltonian and fix a value for the Klein factors, we require that Klein factors 
commute. This is found not be the case here. Since Klein factors satisfy {'7iQ:, Jy^/s} = 25ij8afj, the bosonization of c^^- 
and are independent of the choice of Klein factors, but for i^j, cf^- introduces the Klein factor rjisrjjsrjistljs which 
does not commute with any of the 3-band Klein factors. Furthermore, [r]2sVis'>l2s'']3s,ilis'']3s] 7^ 0. Usually in cases such 
as this one needs to introduce special gauge fields to proceed further with bosonization (if it is possible). Curiously, 
this problem does not seem to arise for ladders with an even number of legs (provided one stays away from dopings 
where 3-leg umklapp scattering is relevant-ie. close to half filling). 

b) N = 4- Bosonization of the relevant 4-band umklapp terms now proceeds without difficulty. Klein factors 



3-leg coupling (in Ui22s) 


4-leg coupling (in C23) 


consistent? 


U1223 = 1 


U?234 = -0.66 
u?324 = -0.69 


yes 


c^j = -0.03 




no 


^2 = 1-0 


f?2 = -0.34 
f?3 = -0.32 


maybe 


f|2 = 0.93 


fj2 = -0.38 
^3 = -0.43 


maybe 


C?2= 1-0 


c?2= 0.002 
c?2= 5x10"'^ 


no 


cf2= 0.91 


CJ2 =0.003 
cl3=-0.0001 


no 


Ui223= -0-26 


Uf234=0.37 
U?324= 0.38 


yes 


0^2= -0.01 


0^3=0.54 


no 


c?i= -0.005 


c?i=-0.11 
c|2= -0.55 


no 


cf3= -0.006 


cf4=-0.18 
ch= -0.59 


no 


u^2i3= -0-24 


Um23= -0.0003 


no 



TABLE V: Correspondence between fixed points: 3-leg ladder RG vs. 4-leg ladder RG with t'=0.1 t. Numbers are relative to 
the respective largest couplings. 



commute, so that one is able to choose the gauge r/is'?2s''?3o-?74CT' = 1 to yield a contribution: 

7A2 E {2<.fe* cos(0(!)+ + (-i)vSt+) cos(<2i, + {-iyei%) 

+ Kjki - cos(</'S+ + (-i)vS'2u) ^os{e<fl,. + {-i)¥flu) 

+ + 4< .,,) cos{cl>'-;i^ + (-l)VS'2i;) cos(eJi\^. + {-l)¥^lu)}, (50) 
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where we have defined 0^^+ = (</)ai +<^a2 + <^a3 + <^a4) and {(j),0)a-ij = y^{{(j),9)ai + {(l),0)aj - {(l),0)a] - {(l>,0)ai)- 
Despite the simphcity of the Klein factors for the 4-leg ladder, it is quite likely that difl^erent physics governs this 
even leg ladder at half-filling than that of the 3-leg ladder, when one considers the RG ratios at which one arrives. 
These are outlined in Table V. One also finds that the 4-band Cooper terms 0^234 and c^234 remain small, so that the 
ground state is not the antiferromagnet. 

c) Large N. In order for the antiferromagnetic terms present at t'=0 and half-filling to grow at large N for t'7^0, a 
large number of umklapp scattering processes need to be relevant. That is, we need to satisfy both the constraints: 
Tlfj^i ^Fj = ^ and the double condition: kpj + kpi + kp-^ + kp^ = 2n, kpk + kpi + kpi + kpj. = 2tt for i ^ j, k. For 
the frustrated Hubbard ladder, no two 2-band couplings are relevant at the same doping. For the splitting of 2-band 
levels presuming no fully filled bands, this double condition (which implies that kpj + kpj = kpk + kp^) seems to now 
require k=j (or at least that they're very close), which is not allowed if it remains a 4-band term. It seems increasingly 
likely that if t' remains finite in the large N limit, we will kill these terms at weak coupling (and the associated AF)! 
Whether such terms may be able to survive for k ~ j at strong U is a matter for future work although it may have 
been partly addressed by work of Honerkamp et alS and Kashima et aliS.. It should be noted that because of the 
fixed nature of the y-component of the momentum in this approach, one ends up treating an asymmetric set of points 
along the effective 2D Fermi surface. When N becomes large this is particularly a problem in that if t' remains fixed, 
as N increases one starts to lose bands to complete filling-that is that the Fermi surface will no longer cross all bands. 

APPENDIX E: BOSONIZATION AFTER BAND 2 MOTT 

Our starting point is 

Ho + Hu„. = (1 _ 5\ (g^^^)2 + !^ ('i _ 5 \ (5^^^)2 ^ Iks^ 

Z \ IT / Z \ IT / TT 

+ — (1 + -j (TTpl) +— (l + -j(^p3) - -^TTplTTpS + (51) 

where TTp is the conjugate field to (ftp as usual, and the spin terms have not been written as the unrotated basis is 
already diagonal for them. 

Replacing the linear combination (tt =f g) "-^^7^-^^ = v^zp, we find that the eigenvectors are given by 




After diagonalization, the kinetic contribution to the Hamiltonian is thus. 



^0 + Hkin 



\ [^-^^ii - f ) ± ^(1 - ^yC-^^Y + i^r] (50,,.)^ 

1 I Vpi+ Vp3 



(1 + J) ^ ^(1 + J).(!^£l_^). + (% ) + 



(53) 



such that the (charge) Luttinger liquid parameters are simple to read off and given in Eq. 24 and 25. 



24 



It would be interesting, then, to bosonizc our order parameter to look for quasi-lD manifestations of phase coherence 
setting in. To do this, it is helpful to re-write our bands at low energy in terms of the linearized spectrum once more, 
and to add the Hermitian conjugate. Then Ai becomes: 

Ai - -^^ +hc - ^^^^'^ ' 



2r/iTr;u(-l)- 
wa 

wa 



Tra 



(cos(V27r(0ai - 0pi)) + cos(V27r((/)ai + Op J)) 



cos(v27r0pi) 



cos 



/ 



v 



/ 



Vd+ 



(54) 



where in the second line we have used the pinning of (j}^^ to set cos(27rn) = 1. Likewise, 



2??3T??3i(-l)' 
7ra 



■ cos 



\ 



w 



Op+b - 



(55) 



The pinning in the charge channel of the operator 9p- then allows us to re-express the arguments of the cosines to 
yield: 



/ 



TTO 



V 



Vd+-1 + \J^ + v; 



,2 

'd+ 



Vd+ 



7p-b — 



2 

d+ 



^p+b 



J 

(56) 

where the forms are identical save for a + sign inside the cosine for the band 1 correlator, and their respective Klein 

factors. Notice that Ai and A3 have been defined above in such a way as to absorb the overall negative sign between 
them expressing their d-wave character. Then the magnitude of the superconducting correlator is: 



2 pK ^ 

< Ai(a;)A3(0) > | = {—f < JJ cos(a7rp + J 



d+ 



-.ep-b{x,Q) 



d+ 



-.ep+b{x,m >, 



(57) 



where in taking the magnitude we can suppress the Klein factors and need not differentiate between sine and cosine 
for p odd and p even for the correlation function. By charge conservation, we know that only terms e*^'"'''e~*^'"''' 
survive, and furthermore, we expect that [Op-b, Op+b] = 0, so that this simplifies to yield: 



aA 



2 )2 
d+' 



2 s2 



and we know that since Kp±i, ^ 1, 9p±h are not canonical fields, so we need to transform to such a basis. We can express 
our results in terms of canonical fields {(fip+b^Op+b) = {\/ Kp±b(j)p±b, —7= — Qp±b) by replacing 6'p±b = y/K^Op±b. 

Simultaneously expanding 6{x) = (j)L{x) — <Pr{x), we arrive at: 



{nay 



< e 



y/Kp-t{4>Lp-b(x)-4>Rf,-b{x)) 



ix/f((^)^(o)) 



25 



"d+- 



y/Kp+b{^Lp+b{x)-^Rp+b(x)) 



X e \V^^'"^+"V^'V e^y¥«-)-(°))+/..c.>, (59) 



and since < e''3'^^-(«)e~*^'^-^(*') > - (4)7^ and < e''5'^«(^)e-*'''^-"(°) > - (i)^, this reduces to Eq. 29 
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